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Abstract 

In this paper we calculate the weak form factors of the decays Bc{B) — > Dlu 
by using the chiral current correlator within the framework of the QCD light-cone 
sum rules (LCSR). The expressions of the form factors only depend on the leading 
twist distribution amplitude (DA) of the D meson. Three models of the L'-meson 
distribution amplitude are employed and the calculated form factor Fb^^d{^) is 
given. Our prediction, by using the D-meson distribution amplitude with the ex- 
ponential suppression at the end points, is compatible with other approaches, and 
favors the three-points sum rules (3PSR) approach with the Coulumb corrections 
included. 
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1 Introduction 



The CDF Collaboration reported on the observation of the bottom-charm Be meson 
at Fermilab [1] in the semileptonic decay mode Be —>■ J/ip + 1 + u with the J/ if) decaying 
into muon pairs in 1998. Values for the mass and the lifetime of the B^ meson were 
given as M{Be) = 6.40 ± 0.39 ± 0.13 GeV and t{Bc) = 0.46l[]i^(stat) ± 0.03(syst) ps. 
Recently, CDF reported first Run II evidence for the Be meson in the fully reconstructed 
decay channel Be J/ip + tx with J/ifj ^ [2\. The mass value quoted for this 

decay channel is 6.2857±0.0053(stat)±0.0012(syst) GeV with errors significantly smaller 
than in the first measurement. Also DO has observed the Be in the semileptonic mode 
Be J/ip + n + X and reported preliminary evidence that M{Be) = 5.95lo;i3 ± 0.34 
GeV and r(5,) = 0.45lg:lg ± 0.12 ps [3]. 

The Be decays, at first, calculated in the potential models (PM) [U [5], wherein the 
variation of techniques results in close estimates after the adjustment on the semileptonic 
decays of B mesons. The Operator Product Expansion (OPE) evaluation of inclusive 
decays gave the lifetime and widths [B] , which agree well with PM, if one sums up the 
dominating exclusive modes. That was quite unexpected, when the sum rules (SR) of 
QCD results in the semileptonic Be widths [7], which are one order of magnitude less 
than those of PM and OPE. The reason may be the valuable role of Coulomb corrections, 
that implies the summation of ag/v corrections significant in the heavy quarkonia, i.e. 
in the Be 0. 

In the recent paper [S], we calculate the form factor for B^Dlu transitions within 
the framework of QCD light-cone sum rules (LCSR). In the velocity transfer region 
1.14 < V ■ v' < 1.59, which renders the OPE near light-cone = go effectively, the 
yielding behavior of form factor is in agreement with that extracted from the data on 
B — »• DW, within the error. In the larger recoil region 1.35 < v-v' < 1.59, the results are 
observed consistent with those of perturbative QCD (pQCD). In this paper we calculate 
the form factor of the semileptonic decay Be^DW, which also depends on the D-meson 
DA. However, due to the different feature of the two process, the c quark is a spectator 
in the decay Be^DW and the c quark comes from the b quark decay in the process 
B DW, these two form factors are sensitive to the shape of the DA in two different 
regions. Combining the information in the two process, we can find which model is 
more suitable for describing the D meson. Similar to the case of B ttIu, the LCSR 
approach for the Be^DW form factor is reliable only in the region < < 15GeV^. 
we extrapolate the result to the whole region and give the decay width and branching 
ratio for the semileptonic decay. 

This paper is organized as follows. In the following section we derive the LCSRs for 
the form factor of Be —>■ DW. A discussion of the DA models for the D meson is given 
in section 3. Section 4 is devoted to the numerical analysis and comparison with other 
approaches. The last section is reserved for summary. 
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2 LCSRs for the B^{B)^D Form Factors 

The Be ^ D weak form factors f{q^) and f{q^) are usually defined as: 

<D{p)\u^^h\B,{p + q)> = 2f{q^)p, + f{q^)q,, (1) 

with q being the momentum transfer. 

To achieve a LCSR estimate of fiq"^), we follow Ref.fTO] and use the following chiral 
current correlator n^(p, g): 

^^.{p,q)=^ J d''xe'P^<D{p)\T{u{x)j,{l + j,)b{x)Xm^ + lMO)}\0> 

= {p + qf)p^ + ti{q\ (p + qf)q^, (2) 

First, we express the hadronic representation for the correlator. This can be done by 
inserting the complete intermediate states with the same quantum numbers as the current 
operator hi{l + 75)0. Isolating the pole contribution due to the lowest pseudoscalar Be 
meson, we have the hadronic representation in the following: 

nf (p, q) = n^(g2, (p + qf)p^ + fi" (p + qf)q^ 
^ <D\u-ff,b\Be><Be\bi-f5c\0> 

^B, -{p + qf 

^ ^ <D\u^,{l + ^,)h\B^><B^\U{l + 75)c|0> 



H 



m 



|« - (p + qy 



Note that the intermediate states B^ contain not only the pseudoscalar resonance of 
masses greater than rriB^, but also the scalar resonances with = 0^, corresponding 
to the operator be. With Eq.([T]) and the definition of the decay constant fs^ of the B^ 
meson 

<5c|6i75c|0> = trIJeJ {rrib + rrie}, (4) 

and expressing the contributions of higher resonances and continuum states in a form of 
dispersion integration, the invariant amplitudes 11^ and 11^ read, 

n^[g', {P + q?] = 7 '^fif^fBjB. ^ r pHs) , ^ subtractions, (5) 



{rrib + rric) {m%^ - {p + qf) J so s - {p + q) 
and 

n^^Ig', {P + q?] = 7 f^'i'^'^BjB P'^is) subtractions, (6) 

^ ^ ^ (mfe + m,)(m|^ - {p + qf) Jso s - {p + qf ' ^ ' 

where the threshold parameter sq should be set near the squared mass of the lowest scalar 
Be meson, the spectral densities p^{s) and p^{s) can be approximated by invoking the 
quark-hadron duality ansatz 

p''is)ip^is)) = p«^^(.)(p«^^(.))^(s - .0). (7) 
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On the other hand, we need to calculate the correlator in QCD theory to obtain the 
desired sum rule result. In fact, there is an effective kinematical region which makes 
OPE applicable: {p + qY — m^-CO for the hd channel and q^<ml — 2KQCD'mb for the 
momentum transfer. For the present purpose, it is sufficient to consider the invariant 
amplitude n(g^, {p+qY) which contains the desired form factor. The leading contribution 
is derived easily by contracting the 6— quark operators to a free propagator: 

< 0|T6(x)6(0)|0 > = / -f^e-^'^^li^. (8) 
J (27r) - ml 

Substituting Eq.® into Eq.([2]), we have the two-particle contribution to the correlator, 

r d^rd'^k 1 
nf ) = -2m,z J ^^e^^^-'^)-^^— -^<D(p)|Tc(x)7,75C^(0)|0 >. (9) 

An important observation, as in Ref. [10]. is that only the leading non-local matrix ele- 
ment <D{p)\u{x)'j^'j5c{0)\0> contributions to the correlator, while the nonlocal matrix 
elements <D{p)\u{x)i'y5c{0)\0> and <D {p)\u{x)a ^^'y5c{0)\0> whose leading terms are of 
twist 3, disappear from the sum rule. Proceeding to Eq.([9]), we can expand the nonlocal 
matrix element <D(p)|T'u(x)7^ 750(0) |0> as 

<D{p)\Tu{x)'j^'y5c{0)\0> = -ip^fo [ (iMe™^^v^£i(M) + higher twist terms, (10) 

where ^d{u) is the twist-2 DA of the D meson with u = 1 — u being the longitudinal 
momentum fraction carried by the c quark, those DA's entering the higher-twist terms 
are of at least twist 4. The use of Eq. (fTO|) yields 

n^^'') (p + qf] = 2fDmb C du „ ^^'^^^ — - + higher twist terms. (11) 

JO mf — [up + qY 

Invoking a correction term due to the interaction of the h quark with a background 
field gluon into Eq. fllll) . the three-particle contribution n^f'^^-* is achievable. However, the 
practical calculation shows that the corresponding matrix element whose leading term is 
of twist 3 also vanishes. Thus, if we work to the twist-3 accuracy, only the leading twist 
DA v^z) is needed to yield a LCSR prediction. 

Furthermore, we carry out the subtraction procedure of the continuum spectrum, 
make the Borel transformations with respect to [p + qY in the hadronic and the QCD 
expressions, and then equate them. Finally, we get the LCSR for /(g^): 



fs ^n{q') = ^^(^^ + ^-)/^e-L/^^^ ^exp 
where 



— (1 — 'u)(g^ — um 



D) 



(12) 



{so" - (Y - m%Y + Am%{ml - q^) - (sf" - - m|,) 

D 



A V ^ uj ' u\ --o ^ I \-u 1 --uj /I Q\ 

« ' 
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and = m\) has been used. 

The LCSR for the form factor of S — Dli) has been derived in Ref. |9|, here we just 
give the result: 



Tb^d{v ■ v') 



2ml 



I mp fp ^rnl/M^ 



where 



A, 



{niB + mp)mB V ttlb Jb 



/■I du 
X — exp 
JAb u 



ml — (1 — u){q'^ — um^)) 



Sq — g2 _ mj^y + Am\,{ml — q^) — {sq — — 



m 



D) 



2m% 



(14) 



(15) 



3 i)-meson Distribution Amplitude 

Now let's do a discussion on an important nonperturbative parameter appearing 
in the LCSRs, the leading twist DA of D-meson, Lpp{x). We reexamine the D-meson 
distribution amplitude since we missed a factor of \/2 for the decay constant fp in 
determining the coefficients of the DA model [9J. 

The D meson is composed of the heavy quark c and the light anti-quark q. The lon- 
gitudinal momentum distribution should be asymmetry and the peak of the distribution 
should be approximately at x ~ m^mp ^ 0.7. According to the definition in Eq. lfTOj) . 
<fp{x) satisfies the normalization condition 



dxipp{x) 



In the pQCD calculations [TT], a simple model (we call model I) is adopted as 



^^ff{x) = 6x(l - x)(l - C,(l - 2x)) 



(16) 



(17) 



which is based on the expansion of the Gegenbauer polynomials. Eq. (ll7l) has a free 
parameter Cd which ranges from to 1, and is supposed to approximate 0.7 in order to 
get consistent results with experiments [H]. Thus we simply take Cd = 0.7. 

On the other hand, it was suggested in that the light-cone wave function of the 
Z^-meson be taken as: 



ipp{x,k±) = Apexp 



-hi 



'k]_ + ml + m; 



+ 



X 



X 



(18) 



which is derived from the Brosky-Huang-Lepage(BHL) prescription [13]. One constraint 
on the wave function is from the leptonic decay process D — > fiu: 



(19) 
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Here the conventional definition of tfie decay constant fo has been used, so Eq. ( fT9i) differs 
from that in Ref . fl2] by a factor of \/2- Another constraint comes from an estimation of 
the probabihty of finding the \qq > Fock state in the D meson: 



D 



dx 



■|^D(x,kx)P 



(20) 



As discussed in Ref.[T2], Pd ^ 0.8 is a good approximation for the D meson. Based 
on these two constraints, the parameters and bj^ can be fixed. Taking Pd ^ 0.8, 
fo = 222.6MeV, = l.SGeV and = 0.35GeV, we have Ad = 225GeV-\ bl = 
0.580GeV-^ ipo {x, k_|_) can be related to the normalized DA foix) by the definition: 



2^6 fd'^k 
7 



D 



'-ipoix, kx). 



Substituting Eq. ffTSj) into Eq. (l2T!) . we have a model of the DA(model II) 



Vqa 



D 



Svr^ fo bl 



■x{l — x) exp 



-hi 



xm?^ + (1 



X 1 



X] 



(21) 



(22) 



Furthermore, as argued in Ref.[14j, a more complete form of the light-cone wave 
function should include the Melosh rotation effect in spin space: 



^D(a;, k_L) = xd{x, ki_)A^jj exp 
with the Melosh factor, 

XD(a;, k_L) = 



'k^ + ml 



D 



X 



+ 



k]j-m|^ 
1 — X 



(1 — x)mc + XTUd 



kj_ + ((1 — x)mc + xmdY 



(23) 



(24) 



It can be seen from Eq. flMj) that XD{.Xj\ii_) — + 1 as rric oo, since there is no spin 
interaction between the two quarks in the heavy-flavor meson, ie., the spin of the heavy 
constituent decouples from the gluon field, in the heavy quark limit [15] . However the 
c-quark is not heavy enough to neglect the Melosh factor. After integration over k_L the 
full form of D meson DA can be achieved (model III): 



XI 




exp 



^/ 2 {xml + {l-x) 
^ x(l 



x] 



(25) 

where y = xrrid -|- (1 — x)mc and the error function Erf{x) is defined as Erf{x) = 
^ exp{—t'^)dt. Using the same constraints as in Eq. ffTOl) and fl20l) . the parameters A{) 

and bi are fixed as A^ = 209GeV"^ and b{j^ = 0.540GeV~l 

In this paper we will employ the above three kinds of models to do numerical calcu- 
lation. All these DA's of the D-meson are plotted in Fig.([l]) for a comparison. It can be 
seen that although they all have a maximum at x ~ 0.65, the shape of them are rather 
different, especial in the region < x < 0.3 and 0.5 < x < 0.8. 
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Table 1: Parameter sets for /^^ and /b, Sq" and Sq for fs^ and /s respectively; rrib, fs^ 
and are given in GeV, s^" and Sq in GeV^. 





mb 


Be 
So 






fs 


set 1 


4.6 


43.0 


0.243 


30.7 


0.145 


set 2 


4.7 


42.0 


0.189 


30.2 


0.117 


set 3 


4.8 


41.2 


0.137 


29.8 


0.090 



4 Numerical Results and Discussion 

Apart from the DA of the D meson, the decay constant of i?c-meson fs^ is among 
the important nonperturbative inputs. For consistency, we use the following corrector 

K{q^) = i y"rfW«^ < 0|c(a;)(l + 75)&(a;),K0)(l -75)c(0)|0 >, (26) 

to recalculate it in the two-point sum rules. The calculation should be limited to leading 
order in QCD, since the QCD radiative corrections to the sum rule for fB^^nio^) are not 
taken into account. The value of the threshold parameter Sq" is determined by a best fit 
requirement in the region 8GeV^<M^<12GeV^, where is the corresponding Borel 
parameter. The same procedure is performed for /s, in almost the same Borel "window" . 
The results are listed in TabiH As we have ignored all the radiation corrections, we don't 
expect our values of fs^ and fs to be good predictions of that quantity. We use the 
same threshold parameters for the corresponding form factors in the LCSRs, except 
for the Borel parameter M^fj, which is taken as ~ M'^/<u>, with <u> been 

the average momentum faraction involved. It turns out that the form factors depend 
little on Ml(j in the region 15 < < 20. The other input parameters are taken as 
TUB = 5.279GeV, mz) = 1.869GeV,mB, = 6.286GeV. 

With these inputs, we can carry out the numerical analysis. In particular, we redo 
the previous calculation for B DW in Ref . [9] and show the corresponding form factor 
Tb^d{v ■ v') in Fig. ([2]). The result for the form factor of DW is given in Fig. ([3]). 

For Tb^d{v ■ v'), similar results can be obtained by applying the various model DA's at 
large recoil region t> ■ t>' ~ 1.59, ie., ~ 0, but rather different values at the zero recoil 
point = q^ax- can be understandable easily from the involved region of the DA. 
While = corresponds to ~ 0.75 according to Eq. (fTB]l . = corresponds 
to ~ 0.6, and the models of the Z)-meson DA in the region 0.5 < x < 0.8 are 
rather different. However, the LCSR result at the zero recoil point (g^ = q^ax) is less 
reliable, we cannot get a final conclusion from the difference of the form factor at this 
point. Fortunately, the case for B^ — »• DW is just opposite, which can be seen from 
Fig-©- There is a big difference of the form factor at the point g^ = 0. A detailed 
comparison for the form factor at this point with other approaches is shown in Tab.©. 
The big difference between model I and others comes from the different contributions of 
the DA's in the involved region < x < 0.45. Due to the exponential suppression at the 
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Table 2: Form factor /^^^/^(O) of DW calculated with different kinds of D-meson 

DA's, in comparison with that of the 3-Points Sum Rule (3PSR) without [7] and with 
[16] the Coulumb corrections and Potential Model (PM) [16] 





medel I 


model II model III 


3PSR [7] 


3PSR [16J 


PM [16J 


/b.^d(O) 


0.55 


0.25 0.28 


0.13 ±0.05 


0.32 


0.29 



Table 3: Form factor fB^^niQ^) in a three-parameter fit (1271) . The three rows correspond 
to the calculated form factors using different sets of parameters, respectively. 





/(O) 


«/ 


bf 


set 1 


0.288 


3.79 


4.23 


set 2 


0.283 


3.92 


4.47 


set 3 


0.288 


4.03 


4.77 



end points, the results from model II and III are much smaller than that from model I, 
and are consistent with the 3PSR results with the Coulumb corrections included, and 
the PM result. It can also be seen from Fig. ([2]) and Fig. ([3]) that, in both cases, model II 
and III actually differ little, which means that the influence of the Melosh factor is not 
so important due to the heavy c quark. 

The calculated form factor for B^^ —>■ Dlv can be fitted excellently in the calculated 
region < < 15GeV^ by the parametrization: 

fBc^oiQ^) = 77 — 2^ } 1 I o I — TT?- i'^'^) 

The values of fBc^oiO), af and 6/ are listed in Table [3l 

Extrapolating the form factor to the whole kinetic region < < {ms^ — ^ 
19.5GeV^ using this parametrization, we get: 

r(5, ^ Dlv) = (0.197 ± 0.013) X lO^^^GeV, (28) 

and 

BR{Bc Dlv) = (1.35 ± 0.05) x 10"^ (29) 

where t{Bc) = 0.45ps and Vub = 0.0037 have been used. The central values are calculated 
by using the parameters set 2, while the upper and lower bounds are given by using 
set 3 and set 1 respectively. Our result for the branching ratio is much larger than 
BR{Bc Dlv) = 0.4 X 10~^ from Ref.[16], they employed a simple pole approximation 
to extrapolate the form factor to the whole region. It is also much larger that the PM 
result BR{B^ Dlv) = 0.35 x 10^^ [17J, and the result BR{B^ Dlv) = 0.6 x 10"^ 
from Ref.dHI. 
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Figure 1: Different kinds of D-meson DA's, solid and dashed curves correspond to model 
III and II, while the dotted line expresses model I. 




Figure 2: J-'b^d as a function of the velocity transfer (with the parameters in the set 2). 
The thin lines expresses the experiment fits results, the solid line represents the central 
values, the dashed(dash-dotted) lines give the bounds from the linear(quadratic) fits. 
The thick fines correspond to our results, with the solid, dashed and dash-dotted fines 
for model III, II and I respectively. 
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5 Summary 



The Be meson has been observed by the CDF and DO groups in the different channels. 
In this paper we study the weak form factor of the decay process Bc{B)—*Dli' by using 
the chiral current correlator within the framework of the QCD light-cone sum rules, 
which is similar to the approach for the weak form factor fBTxiff) in Ref.|10]. The 
calculated form factors depend on the distribution amplitude of the D meson, and we 
employ the three different models for the D meson. It has been shown that the results 
using the model with a exponential suppression at the end points are consistent with 
other approaches. Our results can also confirm the including of the Coulumb corrections 
in the 3PSR calculations for the semileptonic decay B^. DW. In the LCSRs for the 
form factors of Bc{B)^Dli>, the involved region of the D meson distribution amplitude 
is rather different. Combining the information in the two process, we can find which 
model is more suitable for describing the D meson. 

We have made a parametrization fl271) to the form factor by fitting our calculation 
in the region < < 15GeV^ , and the decay width and the branching ratio of 
the process Bc^DW have been calculated. It has been shown that T{Bc —>■ DW) = 
(0.197 ±0.013) X 10-^^GeV and BR{Bc DW) = (1.35 ±0.05) x 10"^ The results are 
different from other approaches. It will be expected to test the different predictions in 
the coming LHC experiments. 



References 

[1] F. Abe et al. (CDF Collaboration), Phys. Rev. D58, 112004 (1998); iMd., Phys. 
Rev. Lett. 81, 2432 (1998). 



[2] D. Acosta et al. (CDF Collaboration), hep-ex/0505076 



[3] M.D. Corcoran (CDF Collaboration), lhep-ex/0506061l 
[4] M. Lusignoh, M. Masetti, Z. Phys. C51, 549 (1991). 

[5] V.V. Kiselev, Mod. Phys. Lett. AlO, 1049 (1995); V.V. Kiselev, Int. J. Mod. Phys. 
A9, 4987 (1994). 

[6] I. Bigi, Phys. Lett. B371, 105 (1996); M.beneke, G.Buchalla, Phys. Rev. D53, 
4991 (1996). 

[7] R Colangelo, G. Nardulh, N. Paver, Z.Phys. C57, 43 (1993); E. Bagan et al, 
Z.Phys. C64, 57 (1994). 

[8] V.V. Kiselev, A.K. Likhoded, A.V. Tkabladze, Phys. Atom. Nucl. 56, 643 (1993), 
Yad. Fiz. 56, 128 (1993); V.V. Kiselev, A.V. Tkabladze, Phys. Rev. D48, 5208 
(1993); V.V. Kiselev, A.K. Likhoded, A.I. Onishchenko, Nucl.Phys. B569, 473 
(2000). 



9 



[9] F. Zuo, Z.H. Li, T. Huang, Phys. Lett. B641, 177 (2006). 



[10] T. Huang, Z. H. Li, and X. Y. Wu, Phys. Rev. D63, 094991 (2001). 
[11] T. Kurimoto, H. N. Li and A. L Sanda, Phys. Rev. D67, 054028 (2003). 
[12] X. H. Guo and T. Huang, Phys. Rev. D43, 2931 (1991). 

[13] S. J. Brodsky, T. Huang and G. P. Lepage, in Particles and Fields-2, Proceedings 
of the Banff Summer Institute, Banff, Alberta, 1981, edited by A. Z. Capri and A. 
N. Kamal (Plenum, New York, 1983), P143; G. P. Lepage, S. J. Brodsky, T. Huang, 
and P. B. Mackenize, ibid., p83; T. Huang, in Proceedings of XXth International 
Conference on High Energy Physics, Madison, Wisconsin, 1980, edited by L. Durand 
and L. G. Pondrom, AIP Conf. Proc. No. 69 (AIP, New York, 1981), plOOO. 

[14] T. Huang, B. Q. Ma, and Q. X. Shen, Phys. Rev. D49, 1490 (1994). 

[15] N. Isgur and M. B. Wise, Phys. Lett. B232, 113 (1989); ibid. 237, 527 (1989). 



[16] V.V. Kiselev, hep-ph/0211021 



[17] M.A.Ivanov, J.G. Korner, P. Santorelli, Phys. Rev. D73, 054024 (2006). 
[18] C.H. Chang and Y.Q. Chen, Phys. Rev. D49, 3399 (1994). 



10 



form factor 
2.5 r 



2.0 - 




0.0 I ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 

0.0 2.5 5.0 7.5 10.0 12.5 15.0 

q'(GeV^) 



Figure 3: /^^^^(g^) calculated by using different kinds of DA models. The solid, dashed 
and dash-dotted lines correspond to Model III, II and I respectively. Here the threshold 
parameter set 2 has been used. 
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